1/1 


^^M5^8^0VMR«I^Rfl^UR^0UQHSESS^^BSHIN^0^HI^EflTTLE 
DEPT  OF  MECHANICAL  ENGINEERING  A  S  KOBAVASHI  ET  AL. 
AUG  85  UHA/DHE/TR-8S/S3  NG8G14-8S-K-ei87 


UNCLASSIFIED 


F/G  11/6 


NL 


line  FILE  COPY  AD-A159  382 


Office  of  Naval  Research 


Contract  N00014-85-K-0187 
Technical  Report  No.  UWA/DME/TR-85/53 
DYNAMIC  FRACTURE  TOUGHNESS 


by 


A.  S.  Kobayashi ,  M.  Ramulu,  M.  S.  Dadkhah,  K.-H.  Yang  and  B.  S.-J.  Kang 


August  1985 


The  research  reported  in  this  technical  report  was  made  possible  through 
support  extended  to  the  Department  of  Mechanical  Engineering,  University 
of  Washington,  by  the  Office  of  Naval  Research  under  Contract  N00014-85-K-0187. 
Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the  United 
States  Government. 


piSnUBimON  STATEMENT 

Apprared  foe  public  nlcoMl  ^.f 
PiittitnittnB  UnJiiaited 


DTIC 


SELECTE 
SEP2  5e85 


B 


Department  of  Mechanical  Engineering 


College  of  Engineering 
University  of  Washington 


Dynamic  Fracture  Toughness 


A.  S.  Kobayashi,  N.  Ramulu,  M.  S.  Dadkhah,  K.-H.  Yang  and  B.  S.-J.  Kang 

University  of  Washington 
Department  of  Mechanical  Engineering 
Seattle,  Washington  98195 


ABSTRACT 

Dynamic  fracture  toughness  versus  crack  velocity  relations  of  Homalite-100, 
polycarbonate,  hardened  4340  steel  and  reaction  bonded  silicon  nitride  are 
reviewed  and  discrepancies  with  published  data  and  their  probable  causes  are 
discussed.  Data  scatter  in  published  data  are  attributed  in  part  to  the 
observed  fluctuations  in  crack  velocities.  The  results  reaffirmed  our  previous 
conclusion  that  the  dynamic  fracture  toughness  versus  crack  velocity  relation 
is  specimen  dependent  and  that  the  dynamic  crack  arrest  stress  intensity  factor 
is  not  a  unique  material  property. 

INTRODUCTION 

Since  Wells  and  Post  fll,  with  the  help  of  Irwin  121,  determined  the  crack 
driving  force,  i.e.  the  dynamic  stress  Intensity  factor,  and  the  crack  velocity 
in  fracturing  photoelastic  plates,  numerous  attempts  have  been  made  to  relate 
these  two  quantities.  The  dynamic  fracture  community's  interest  in  this 
relation  is  demonstrated  by  the  fact  that  six  out  of  the  seven  review  papers 
dealing  with  the  experimental  aspects  of  dynamic  fracture  mechanics  in  the 
recent  issue  of  the  International  Journal  of  Fracture  [31  refer  to  the 
uniqueness  or  lack  thereof  in  the  dynamic  stress  intensity  factor  versus  crack 
velocity  relation  and/or  in  the  dynamic  crack  arrest  stress  intensity  factor. 

The  survey  paper  by  Dally  et  al  r43  describes  the  major  findings  to  date  and 


indicates  possible  sources  of  experimental  errors  idiich  ray  have  lead  to  the 
current  controversies  on  this  subject. 

The  purpose  of  this  paper  is  to  present  additional  experimental  results, 
some  of  which  were  obtained  by  the  authors  and  their  colleagues  over  the  past 
decade,  on  dynamic  stress  intensity  factor  versus  crack  velocity  relations  in 
the  context  of  the  current  controversy.  Throughout  this  paper,  the  measured/ 
computed  dynamic  stress  intensity  factors  are  referred  to  as  the  dynamic 
fracture  toughness.  Thus  the  driving  force,  i.e.  the  dynamic  stress  intensity 
factor.  Is  tacitly  equated  to  the  material  resistance  to  dynamic  crack  growth, 
i.e.  dynamic  fracture  toughness. 

DYNAMIC  PHOTOELASTIC  RESULTS 

Although  photoelastic  polymers,  such  as  Homalite-100  and  epoxy,  are  not 
primary  structural  material,  dynamic  photoelasticity  and  caustics  have  been 
used  in  the  past  decade  and  half  to  uncover  the  basic  principles  which  govern 
dynamic  fracture  mechanics.  The  dynamic  fracture  toughness,  Kjq,  versus  crack 
velocity,  a,  relation.*?,  which  have  been  obtained  through  extensive  fracture 
testing  of  polymers,  showed  that  the  the  terminal  crack  velocity  is  test 
specimen  dependent  while  the  "near  vertical  stem"  of  these  relations  is  either 
a  unique  C51  or  a  nonunique  C6,71  material  property.  The  latter  is  in 
agreement  with  the  conclusion  derived  by  one  of  the  authors  several  years  ago 
C8.9J.  The  dynamic  photoelastic  data  used  to  support  this  conclusion  has  been 
reevaluated  in  this  paper  by  an  updated  data  processing  procedure  which 
incorporates  higher  order  terras  of  the  dynamic  crack  tip  stress  field. 

Figures  1  and  2  .•?how  the  Kjp  versus  a  relations  for  Homalite-lOO  and 
polycarbonate  fracture  specimens.  No  attempt  was  made  to  fit  an  average  Kjo 
versus  a  curve  through  the  wide  scatter  of  data  generated  from  various  batches 
of  Horaalite-100  and  polycarbonate  sheets  tested  over  a  period  of  ten  years. 


Figure  1  shows  that  the  scatter  bands  about  the  imagined  vertical  stems  of  the 
dynamic  tear  test  (DTT),  single  edge  notched  (SEN),  modified  compact  (M-CT)  and 
wedge-loaded  rectangular  double  cantilever  beam  (HL-RD(ZB)  Homalite-100 
specimens  are  similar  to  those  shown  in  C7].  However,  differences  in  the 
minimum  dynamic  fracture  toughness,  of  the  vertical  stems  of  the  DTT  and 
SEN  specimens  are  larger  than  that  reported  in  C6D.  The  difference  in  Kju,  for 
the  more  ductile  WL-RDCB  and  DTT  polycarbonate  specimens  is  about  10  percent 
and  is  in  agreeement  with  the  general  observation  by  Rosaliis  et  al  ClOl. 

Figure  3  shows  the  Kjq  versus  crack  extension  relations  of  four  SEN 

specimens  subjected  to  different  fixed  grip  loading  condition  £91.  Also  shown 

are  the  corresponding  static  stress  Intensity  factor.  This  figure,  which  is 

similar  to  the  well-publicized  results  of  Kalthoff  et  al  CllD,  demonstrates 

that  the  dynamic  crack  arrest  stress  intensity  factor,  is  a  constant  for 
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the  same  specimen  while  the  static  crack  arrest  Intensity  factor,  Kjj,  ,  varies 
with  the  crack  initiation  condition. 


SCATTER  IN  Kjp  VERSUS  a  RELATION 

Since  the  above  photoelastic  results  are  in  general  agreement  with  the 
caustic  results,  the  published  discrepancies  in  the  Kjq  versus  a  results  cannot 
be  attributed  to  the  differences  in  the  experimental  procedures  alone. 

However,  the  discrepancies  could  be  attributed  in  part  to  the  size  of  the  crack 
tip  region  used  for  data  reduction  in  the  presence  of  stress  wave  effects  C123. 
The  caustic  method  by  definition  and  the  authors'  photoelastic  method  by  choice 
had  restricted  the  crack  tip  region  to  within  5  ram  of  the  crack  tip  but  outside 
of  the  nonlinear  region  of  about  1  ram  C131  surrounding  the  crack  tip.  The 
dynamic  photoelastic  results  in  C43  are  derived  from  larger  crack  tip  regions 
with  the  use  of  larger  number  of  higher  order  terms  in  the  crack  tip  stress 
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field.  Such  data  reduction  procedure  Hill  yield  accurate  stress  intensity 
factors  under  static  loading.  On  the  otherhand,  the  dynanic  isochroaatics  in  a 
larger  crack  tip  region  Hould  be  less  sensitive  to  snail  perturbations  in  the 
dynamic  stress  intensity  factor  as  shoHn  in  a  previous  numerical  experiment 
CSl.  The  combined  effect  of  the  large  crack  tip  region,  in  Hhich  measurements 
were  made,  and  the  large  fracture  specimens  C41,  which  are  shown  in  the  right 
half  of  Figure  4,  would  minimize  any  oscillations  in  the  Kjq.  In  contrast,  the 
stress  wave  effect  is  more  severe  in  the  smaller  fracture  specimens,  which  is 
shown  in  the  left  half  of  Figure  4,  and  the  resultant  oscillations  in  Kjo  is 
more  readily  detected  when  a  smaller  crack  tip  region  is  used  in  data 
reduction. 

The  experimental  errors  involved  in  crack  velocity  measurements  have  been 
discussed  in  C4,12,153  with  C41  suggesting  the  use  of  ultrasonic  fractography 
C161  for  increased  accuracy.  Such  crack  velocity  measurements  C173  were  made 
on  CT,  SEN,  3-polnt  bend  and  Charpy  polymethyl  methacrylate  (PfflA)  specimens 
approximately  one  half  the  size  of  the  smallest  WL-RDCB  specimen  in  Figure  4. 
The  qualitative  changes  in  the  crack  velocities  with  crack  extension  in  the 
SEN,  CT  and  3-point  bend  specimens  are  similar  to  those  reported  in  [18,19,203, 
respectively.  Moreover,  the  crack  velocities,  which  were  determined  from  the 
discrete  Cranz-Schardin  photographs,  ultrasonic  fractography  and  streaking 
photography  did  not  exhibit  any  unusual  perturbation  in  the  otherwise  gradually 
varying  crack  velocities  in  these  polymeric  materials.  Figure  5  shows  the 
experimental  setup  and  a  typical  streaking  photograph  C193  used  to  determine 
the  continuous  change  in  crack  velocity  in  a  fracturing  polycarbonate  modified 
compact  (M-CT)  specimen.  While  the  crack  velocity  measurements,  which  were 
made  directly  from  the  Cranz-Schardin  photographs  may  not  be  accurate,  the 
results  appear  to  be  in  qualitative  agreements  with  those  obtained  by  the  more 
accurate  ultrasonic  fractography  C173  and  streaking  photography  [193. 


The  small  but  sharp  changes  in  crack  velocities,  which  are  comparable  to 
those  reported  in  tl2D,  were  observed  in  the  Charpy  specimens  C17D  which  were 
subjected  to  severe  stress  wave  effects.  As  will  be  shown  later,  such 
discontinuous  crack  velocities  was  also  observed  in  small  hardened  4340  steel 
and  ceramic  specimens  where  the  stress  wave  effect  is  pronounced. 

Experimental-Numerical  Procedure 

The  crack  tip  state  of  stress  of  a  propagating  crack  in  opaque  or  optically 
insensitive  material  has  been  determined  by  photoelastic  coating  method  C213 
and  the  more  popular  caustic  method.  An  alternate  procedure  is  to  combine 
experimental  and  numerical  techniques  by  using  measured  crack  extension  history 
interactively  with  a  dynamic  finite  difference  or  finite  element  program  in  its 
generation  or  propagation  phase  C223.  The  latter  propagation  analysis  was  used 
extensively  by  the  Battelle  group  C233  to  study  the  dynamic  crack  toughness  and 
arrest  charateristics  of  steel  C243  and  by  one  of  the  authors  and  his  colleague 
to  study  the  nonlinear  fracture  response  of  concrete  C253.  The  former 
generation  analysis  has  been  used  to  study  the  dynamic  fracture  respon.se  of 
glass  C26]  and  reaction  bonded  silicon  nitride  r253. 

The  above  hybrid  experimental-numerical  procedure  was  used  to  determine 
the  Kjq  versus  a  relation  for  4340  steel  hardened  to  Rockwell  C  44.  The 
dynamic  crack  extension  histories  in  four  wedge-loaded  modified  double 
cantilever  beam  specimens  (WL-MDCB),  shown  in  Figure  6  (a),  and  with  a  chevron 
starter  notch  were  measured  by  a  ICRAK-GAGE*  and  FRACTOMAT.*  Figure  7  shows 
typical  crack  extension  records  of  two  fracturing  4340  WL-MDCB  specimens.  The 
initial  and  slower  crack  propagation  in  the  chevron  notch  specimens  is  followed 
by  rapid  crack  propagation  and  subsequent  deceleration.  The  latter  crack 
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deceleration  is  interrupted  by  a  number  of  short  intervals  of  crack  arrest 
where  the  average  time  between  each  crack  arrest  coincides  with  the  average 
transit  time  of  shear  wave  from  the  crack  tip  to  the  lateral  edge  of  the 
specimen  and  back. 

Such  intermittent  crack  propagation  is  more  pronounced  in  the  blunt  notch 
4340  WL-MTDCB  specimen,  which  is  heat  treated  to  a  hardness  of  Rockwell  C  52. 
Figure  8  shows  the  crack  extension  history  with  crack  arrest  Intervals 
indicated  by  arrow  marks.  Such  intermittent  crack  arrests,  as  long  as  20 
microseconds,  were  reported  by  Van  Elst  C28D  and  de  Graaf  C293,  who  used 
streaking  photography  to  record  continuous  crack  extension  in  Robertson  type 
low-carbon  steel  specimens.  Ravi-Chandar  et  al  C73  and  Rosakis  et  al  CIO!  also 
reported  the  presence  of  discontinuous  crack  velocities  in  their  highly 
dynamically  loaded  specimens. 

Returning  to  the  hybrid  experimental-numerical  procedure,  an  average  of 
the  measured  crack  extension  histories,  which  are  shown  in  Figure  9,  without 
crack  arrest  of  four  4340  steel  WL-MTDCB  specimens  was  then  used  to  drive  a 
dynamic  finite  element  code  in  its  generation  mode  and  the  dynamic  fracture 
parameters  were  determined. 

Figure  10  shows  the  Kjq  versus  crack  extension  relation  as  well  as  the 
corresponding  static  stress  intensity  factor  in  this  high  strength  4340  steel 
WL-MTDCB  specimen.  Figure  11  shows  the  Kjq  versus  a  relation  for  this  study 
as  well  as  that  of  Rosakis  et  al  nOl.  The  remarkable  agreement  between  the 
two  independent  results  could  be  due  in  part  to  the  similarities  in  specimen 
geometries. 

Despite  the  differences  in  fCjp  versus  a  relations,  a  vertical  stem  in 
the  Kjq  versus  a  relation  always  existed  in  the  photoelastic  polymers  and  4340 
steel  specimens  discussed  so  far.  However,  limited  dynamic  fracture  studies  of 


extremely  brittle  materials,  such  as  glass  and  structural  ceramics  C26,27,303, 
show  that  ICj^  and  hence  the  vertical  stem  in  the  fCjp  versus  a  curve  does  not 
exits  in  some  materials.  Figure  12  shows  the  Kjq  versus  A  relation  of  reaction 
bonded  silicon  nitride  WL-MTDCB  specimens  loaded  to  fracture  under  both  static 
and  dynamic  conditions.  The  specimen  geometry  is  identical  to  that  shown  in 
Figure  6  (a)  with  Figure  6  (b)  showing  the  dynamic  loading  arrangement.  While 
the  crack  propagating  under  static  loading  had  attempted  to  arrest,  as  shown  in 
Figure  12,  the  same  crack  propagating  under  dynamic  loading  showed  little 
tendency  for  arresting. 

CONCLUSIONS 

As  profoundly  stated  by  many  authors  in  C3D,  the  controversy  regarding  the 
uniqueness  or  lack  thereof  in  the  Kjq  versus  a  relation  is  far  from  being 
settled.  While  available  experimental  results  indicate  that  in  the  absence  of 
stress  wave  effects,  such  as  in  infinitely  large  fracture  specimen  under  benign 
loading,  Kjq  versus  a  relation  may  possess  a  unique  or  a  vertical  stem. 

Such  unique  vertical  stem  is  not  observed  in  dynamic  fracture  specimens  of 
smaller  size  and/or  under  dynamic  loading. 

Comparative  study  of  various  experimental  data  shows  that  the  consistency 
in  data  scatter  cannot  be  totally  attributed  to  experimental  errors  and  that 
the  intermittent  crack  arrest  and  the  discrete  changes  in  crack  velocity  are 
caused  by  the  reflected  stress  wave. 

DISCUSSION 

In  the  pursuit  of  the  above  uniqueness  controversy,  we  pose  the  question 
"for  what  reason?"  The  end  use  of  the  sought  Kjo  versus  a  relation  is  as  the 
fourth  constitutive  equation  for  estimating  the  dynamic  fracture  response  of  an 
elastic  solid.  Limited  numerical  experiments  show  that  the  arrest  crack  length 


of  a  propagating  crack  is  obviuosly  governed  by  Kj,^  C31-343.  For  a  dynamically 
loaded  specimen  or  in  the  presence  of  severe  stress  wave  effects,  however, 
small  differences  in  may  not  cause  large  differences  in  the  arrest  crack 
length  while  the  same  difference  in  Kj^,  may  cause  large  differences  in  arrest 
crack  length  in  the  absence  of  stress  wave  effects. 
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